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We present a generalized dynamical mean-field approach for the nonequilibrium physics of a strongly 
correlated system in the presence of a time-dependent external field. The Keldysh Green's function 
formalism is used to study the nonequilibrium problem. We derive a closed set of self-consistency 
equations in the case of a driving field with frequency fl and wave vector q. We present numerical 
results for the local frequency-dependent Green's function and the self-energy for different values 
of the field amplitude in the case of a uniform external field using the iterated perturbation theory. 
In addition, an expression for the frequency-dependent optical conductivity of the Hubbard model 
with a driving external field is derived. 

PACS numbers: 71.10.Fd, 71.27. +a 
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The study of nonequilibrium many-body phenomena 
is a field of great and still growing interest in today's 
condensed matter physics. Recent pump-probe experi- 
ments on transition metal oxides show unusual effects, 
e.g. a dramatic change in the optical transmission!!!, 
which require a theoretical description of driven quan- 
tum systems. A basic model for theoretical investigations 
of strongly correlated fermion systems is the Hubbard 
model. Within—the scope of the dynamical mean-held 
theory (DMFT)B, a very detailed analysis of the spectral 
properties of the Hubbard model became possible. The 
DMFT approach takes into account local correlations as 
well as high energy excitations and provides a useful tool 
for calculating linear response functions of strongly cor- 
related fermion systemscm. 

In this paper, we present a generalization of the DMFT 
using the Keldysh formalism to study a Hubbard model 
under the influence of a periodic external perturbation. 
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where cJ CT (c i(T ) creates (annihilates) a fermion with spin 
a at site i and = c\ a c i(y is the local density operator. 
The bare half-filled single-band Hubbard model describes 
nearest-neighbor hoppings of electrons on a lattice, inter- 
acting with each other through the on-site Coulomb inter- 
action. The additional explicitly time-dependent term in 
the Hamiltonian describes the external driving field with 
frequency f2, wave vector q = (g, q) and arbitrary field 
amplitude A. 



Vi(t) = A cos (fit - qR,) 
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In particular we are interested in the small-g limit. This 
Hamiltonian describes photodoping of a strongly corre- 
lated material, since the photon wavelength is large com- 
pared to the lattice spacing. For a uniform external field 



(q — 0) the time-dependent perturbation term commutes 
with the Hubbard Hamiltonian, and this problem can be 
solved quite easily, as we will see later. 

For large commensurate wave vectors, e.g. q = ir, 
the d-dimensional lattice can be divided into cluster of 
N = 2n/q sites. The local properties depend on the 
chosen site R of this cluster. Performing the limit of a 
small wave vector is numerically limited by the increas- 
ing number of sites in the cluster. A theory which syste- 
matically includes nonlocal correlations to the DMFT is 
the dynamical cluster approximation (DCA)Ij, where the 
original model is mapped onto a self-consistent embedded 
cluster instead of a single impurity. 

The main technique for the treatment of nonequilihrl- 
um effects is the Keldysh Green's function formalismH'LI. 
The usual time-ordering is replaced by a contour-ordering 
along the Keldysh contour, which follows the real-time 
axis from — oo to +oo and back to — oo. The resulting 
Green's functions matrix G"' 3 , where a,f3 G {±1} speci- 
fies which time argument lies on the upper (+), respec- 
tively lower (— ) branch of the Keldysh contour, defines 
three linearly independent components which provide in- 
formation about the quantum states of the system (re- 
tarded/advanced Green's function) as well as about their 
occupation out of equilibrium within a density matrix 
formulation (Keldysh component). 

The external field is periodic in time and in each spa- 
tial coordinate. Thus, according to a generalized Flo- 
quet theoremH, energy and momentum are conserved up 
to multiples of the external frequency and wave vector, 
respectively. We define a double Fourier transformation, 
expanding the Green's function in its Floquet modes 
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In the noninteracting case (U = 0) and for small q, it is 
possible to calculate the Green's function analytically by 
solving the equation of motion. The weighting coefficient 
of the n-th energy sideband turns out to be the Bessel 
function J n (A/il q (k)) with r2 g (k) = Q — qvt- Thus, 
the noninteracting retarded/advanced Green's function 
is given by 



Jm (-Ag (k) ) J m — n (A q (k) ) 

lu — — mf2 g (k) ± i0+ 



(4) 



where A g (k) = A/0„(k). The noninteracting Keldysh 
component reads 
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(k,w) =-27ri^(l-2/(e k _ mq )) 
xJ m (Ag(k))J m _ n (A g (k)) 



fe(*a, ig)) depending on the chosen site R of the cluster. 
The self-consistency condition is obtained by integrating 
out all degrees of freedom except those of the selected 
site R. Assuming a semicircular bare density of states 
(Bethe lattice) one obtains the self-consistency condition 
for the Fourier transformed Weiss field, which is a matrix 
equation in energy and momentum sidebands and also in 
Keldysh space 



(u> — nf2) S„ 



"2 {^nm+l + $nm-l) 



-a/3t 2 cos((n-m) g )G^„>). (7) 



x5 {ijj — ek - mQ q (k)) . (5) 

The diagrammatic expansion of the Keldysh Green's 
function matrix in the Coulgmb parameter U is mainly 
identical to Feynman theoryQ. The only difference arises 
from the summation over the two time branches. 

We solve the driven Hubbard model using the dyna- 
mical mean-field approximation. Within DMFT the ori- 
ginal lattice problem is mapped onto an effective impu- 
rity problem, which can be solved numerically, using e.g. 
quantum-Monte-Gadp simulations or iterated perturba- 
tion theory (IPT)ISEl The local dynamics at a specific 
site R can be regarded as the interaction of the degrees 
of freedom at this site with an external bath created by 
all other degrees of freedom on other sites. Due to the 
spatial fluctuations for nonzero external wave vector the 
local properties depend on the chosen site R of the clus- 
ter. 

The key approximation of the DMFT is the assump- 
tion of a local self-energy, which-caii be justified in the 
limit of large lattice connectivitjO'E3. For the bare Hub- 
bard model momentum conservation can be disregarded 
at each internal interaction vertex in the limit of infi- 
nite connectivity. An analogous argument for the driven 
Hubbard model yields that at internal vertices momen- 
tum conservation is automatically fulfilled up to multi- 
ples of q. Thus, each internal propagator G"f„ 2 (k, w) — 
^k-niq k-n 2 q( tJ — Tiiil, u> — r^fi) can be replaced by 
Gq%(uj) = J2k ^9,0n( k ) w ) which is a "staggered" average 
over the lattice sites R in the cluster: ^ R e _mqR GR,(w). 
The resulting Dyson equation reads 

G qfin M = 9qfin M 

+ E^m^Lhg^m- (6) 

The mean-field description of the driven Hubbard 
model is determined by the effective Weiss field 



We used the iterated perturbation theory (IPT)I§£2I to 
solve the mean-field equations numerically for the half- 
filled driven Hubbard model. The basic idea of the IPT 
method is to perform a skeleton expansion of the self- 
energy up to second order in the interaction U : 
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x G J q ?o, n+m+ i ("1 + "a) Gt,-m - ^) ■ (8) 



The matrix G g , nm is then obtained from the Dyson equa- 
tion ^) with the noninteracting g q . nm replaced by the 
Weiss field Gq.nm- The self-consistency loop is closed by 
calculating a new Weiss field from equation (0). 

In our numerical calculation, we consider the special 
case of a uniform external field (q = 0), where the time- 
dependent perturbation term commutes with the bare 
Hubbard Hamiltonian and therefore the weighting coeffi- 
cients of the energy sidebands are Bessel functions in the 
interacting case, too. 

For A = 1 and A = 2 (we set f2 = 1), we restricted 
our considerations to the main energy band and the first 
three sidebands, i.e. we calculated the Fourier compo- 
nents Go n (uj) with n = —3, .., 3 (see Fig. Q). For A = 3, 
we considered five sidebands: n — —5, ..,5. These limi- 
tations of the number of sidebands are reasonable, since 
for a given A the Bessel functions J n (A) fall off very fast 
with increasing n. This property strongly depends on the 
magnitude of the field: we have to take into account an 
increasing number of sidebands while enhancing the field 
amplitude. 
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FIG. 1. ImGo^(u>) for U = 3, A = 1, three sidebands. 



FIG. 3. ImE^Lj) for U = 3 and A = 0, 1,2, obtained by 
IPT. 



As an example, the Fourier components of the spectral 
density for A = 1 are given in Fig. [l| From an inverse 
Fourier transformation of the advanced Green's function, 
we obtained that the time-dependent advanced function 
satisfies the usual causality condition. 

In the presence of the uniform external field, we define 
an effective spectral density, given by 



p«*( e ) = IlmGgoH = V Jl(A)-lmG(oj - 1(1), (9) 

I 



where G is the local Green's function for A = 0. Fig. g 
shows how the effective spectral density depends on the 
field amplitude. In comparison with the usual three peak 
structure in the case of no external field (lower and upper 
Hubbard band and the quasiparticle peak), this function 
exhibits a more complex structure for A due to the 
sidebands, weighted by the Bessel functions. The num- 
ber of occuring sidebands depends strongly on A: with 
increasing A, more sidebands become significant in the 
effective spectral density. An additional phenomenon is 
the strong suppression of the main quasiparticle peak for 
certain field amplitudes. 

The effect of the external field on the imaginary part of 
the advanced self-energy on the field amplitude is shown 
m Fig. |. For A ^ 0, we observe the sidebands of the 
typical two-peak structure in the equilibrium case. 




FIG. 2. ImGooH for U = 3 and A = 0, 1,2, obtained by 
IPT. 



Finally, we consider the Keldysh component of the 
Green's function, which contains information about the 
occupation of the states. From Goo ld we define an effec- 
tive distribution function: 



(10) 



For zero temperature and A = 0, / off is the usual Heavi- 
side function. The effective distribution functions for 
A = 1,2,3 are shown in Fig. ^. The modifications due 
to the external field at the energy sidebands of the Fermi 
edge correspond to emission and absorption processes of 
energy quanta f2. Experiments to measure the energy 
distribution function in a stationary out-of-equilibrium 
situation in metallic wires connected to two reservoir elec- 
trodes with a*| applied bias voltage have been done by 
Pothier et al.l 




FIG. 4. Effective distribution function for A = 0, 1, 2, 3 at 
T = 0. 
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The time-dependent DMFT also provides a useful 
method for calculating the optical conductivity of the 
driven Hubbard model. Consider electrons described by 
the driven Hubbard model, which couple to an additional 
external classical vector potential with wave vector p. 
Linear-response theory for noncquilibrium systems can 
be formulated quite simple in terms of Keldysh propaga- 
tors. The paramagnetic current response in e— direction 
is given by 

( J?-™ (p, t)> = f £ v, (k + p, k) catgut, t) (ii) 



with the vertex factor 



. (k + p, k) = e j(k+p)s - 



-zke 



(12) 



Just as in the case of the bare Hubbard Hamiltonian, 
all vertex corrections in the ladder decomposition of the 
current response function vanish in the limit of infinite 
dimensionstl "Ho. Thus, the optical conductivity is deter- 
mined by the elementary particle-hole bubble as in the 
equilibrium case. Assuming isotropic spatial conditions, 
we obtain for the paramagnetic contribution 
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(w' + w) (13) 



where G^ denotes the zeroth order contribution in the 
classical vector potential, T\ is the first Pauli matrix and 
v k = Vkfik- Using the spectral representation of the 
Green's functions, the optical conductivity can be ex- 
pressed in terms of the Fourier components of the spec- 
tral density p l 0n (e) = J t (£) J ; _ n (^)ImG av (e - IQ)/tt and 
the free density of states D(e) (the diamagnetic term can- 
cels the 1 /w-divergence of the paramagnetic part in the 
optical conductivity). 
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J de J du'D(e) 



X (/(u/-/ 2 fi)-/(o/ + cJ-Zifi)) (14) 

This final result for cr(u>) significantly shows the appea- 
rance of the energy sidebands due to the external driving 
field, and for zero driving field it reduces to the well- 
known expression for the optical conductivity of the bare 
Hubbard model. 

We have generalized the DMFT equations using the 
Keldysh formalism to study a Hubbard model under the 
influence of a strong periodic external perturbation with 
arbitrary frequency and amplitude. We found that the 
Green's functions were matrices in the energy and mo- 
mentum sidebands due to quasi-energy and momentum 



conservation modulo the external frequency respectively 
wave vector. 

We have derived a closed set of self-consistency equa- 
tions for the driven Hubbard model wave vector nume- 
rically in the special case of a uniform field using IPT. 
We obtained an effective spectral density and distribu- 
tion function for different values of the field amplitude. 
In addition, we have derived the frequency-dependent op- 
tical conductivity of the driven Hubbard model. 

Solving the problem for a nonzero commensurate q 
would require an even larger numerical effort. One has 
to consider cluster of N = 2ir/q sites and perform time- 
dependent DMFT for all sites R of this cluster. This is 
basically the. same as using dynamical cluster approxima- 
tion (DCA)til in combination with the Keldysh formalism. 
The original driven Hubbard model is mapped onto a self- 
consistent embedded cluster given by the periodic form 
of our Hamiltonian instead of a single impurity. Decreas- 
ing the wave vector requires solving cluster problems of 
increasing size, which is numerically limited because of 
an exponentially growing Hilbert space. 

We are currently working on an application of the time- 
dependent DMFT to describe the local nonequilibrium 
physics of excitonic insulators and Mott insulators un- 
der laser excitation. This method shows great promise 
for being suitable to explain recent pump-probe experi- 
ments on transition metal oxides where unusual effects, 
e.g. a dramatic change in the optical transmission^, have 
been observed. The real-time development of the DOS 
and the gap after the impact of a short intense laser pulse 
is under investigation. 



1 T. Ogasawara et al. Phys. Rev. Lett. 85, 2204 (2000) 

2 A. Georges, G. Kotliar,, W. Krauth, M. J. Rozenberg, Rev. 
of Mod. Phys. 68, No 1, 13 (1996) 

3 Th. Pruschke, D. L. Cox, M. Jarrell, Europhys. Lett. 21, 
No 5, 593 (1993) 

4 Th. Pruschke, D. L. Cox, M. Jarrell, Phys. Rev. B 47, No 
7, 3553 (1993) 

5 M. H. Hettler et al. Phys. Rev. B 58, No. 12 (1998) 

6 L. V. Keldysh Sov. Phys. JETP 20, No 4, 1018 (1965) 

7 J. Rammer, H. Smith, Rev. of Mod. Phys. 58, No 2, 323 
(1986) 

8 P. Hanggi, Quantum transport and dissipation, Driven 
Quantum Systems, (1998) 

9 A. Georges, G. Kotliar Phys. Rev. B 45, No 12, 6479, 
(1992) 

10 E. Muller-Hartmann, Z. Phys. B 76, 211 (1989) 

11 W. Metzner, D. Vollhardt, Phys. Rev. Lett. 62, No 3, 324 
(1989) 

12 E. Muller-Hartmann, Z. Phys. B 74, 506 (1989) 

13 H. Pothier et al. Phys. Rev. Lett. 79, No. 18 (1997) 

14 A. Khurana, Phys. Rev. Lett. 64, No. 16, 1990 (1990) 



4 



